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j^ \ Kinesins are processive motor proteins that move along microtubules in a step- 

wise manner, and their motion is powered by the hydrolysis of ATP. Recent experi- 
P3 ■ ments have investigated the coupling between the individual steps of single kinesin 

o, 

^O^, molecules and ATP hydrolysis, taking explicitly into account forward steps, back- 

ward steps and detachments. A theoretical study of mechanochemical coupling in 

> 

00 ' kinesins, which extends the approach used successfully to describe the dynamics of 

lO . conventional motor proteins, is presented. The possibility of irreversible detach- 

^, 

^^ \ ments of kinesins from the microtubules is also explicitly taken into account. Using 

■^ ' 

>-^ ' the method of first- passage times, experimental data on the mechanochemical cou- 

pling in kinesins are fully described using the simplest two-state model. It is shown 
_^ , that the dwell times for the kinesin to move one step forward or backward, or to 

Q I dissociate irreversibly are the same, although the probabilities of these events are 

O ■ 

different. It is concluded that the current theoretical view, that only the forward 

k> , motion of the motor protein molecule is coupled to ATP hydrolysis, is consistent 

H ; 

C^ \ with all available experimental observations for kinesins. 
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INTRODUCTION 

There are several classes of enzymes, called molecular motor proteins, that are critical 
for many biological processes, but especially they are important for cellular transport and 
motility, cell division, and transfer of genetic information (Lodish et al., 1995; Bray, 2001; 
Howard, 2001). The motor proteins, such as kinesins, myosins, DNA and RNA polymerases, 
move in a stepwise motion along rigid molecular tracks (microtubules, actin filaments, and 
DNA molecules) . The motion of motor proteins is fueled by the hydrolysis of ATP or related 
compounds. However, the exact mechanism of the coupling between the chemical energy of 
hydrolysis and the mechanical motion of motor proteins is still unknown, and it remains one 
of the most important problems in biology. 

Kinesins provide the most convenient system to investigate the mechanochemical coupling 
in motor proteins since biophysical, chemical and mechanical properties of these molecules 
are now well studied at single-molecule level (Howard et al., 1989; Svoboda et al., 1994; 
Schnitzer and Block, 1997; Kojima et al., 1997; Visscher et al., 1999; Schnitzer et al., 2000; 
Nishiyama et al., 2003; Asbury et al., 2003, Yildiz et al., 2003). Conventional kinesins are 
dimeric two-headed molecules, which hydrolyze ATP and move stochastically in 8.2-nm steps 
along the microtubules. These motor proteins can make hundreds of steps before dissociating 
from the microtubules and they can be processive even against the opposing load as high as 
7-8 pN (Visscher et al., 1999; Schnitzer et al., 2000; Nishiyama et al., 2002). Kinesin moves 
preferentially in the forward direction (plus end of microtubules), however, at high loads the 
frequency of backward steps (in the direction of minus end of the microtubule) is increasing 
(Nishiyama et al., 2002). 

In order to understand how the motor proteins function, it is important to investigate 
how the chemical energy of ATP hydrolysis is transformed into the mechanical motion of 
proteins. To approach this fundamental problem, first, several critical questions should be 
answered: 1) How many ATP molecules consumed for each kinesin's step? 2) Are ATP 
molecules hydrolyzed for any step, forward or backward? 3) Is there a futile hydrolysis in 
kinesin motion, i.e., ATP consumption without actual moving of the motor protein? 

In recent experiments (Nishiyama et al., 2002), the mechanism of mechanochemical cou- 
pling in motor proteins has been studied by correlating the forward and backward move- 
ments of single kinesin molecules to the hydrolysis of ATP. Using optical trapping nanometry 



system, the time trajectories of single kinesin molecules have been measured for different 
external forces and for different ATP concentrations. It was found that the dwell times 
before the forward and backward steps are the same at all external forces and at all ATP 
concentrations. A biased Brownian motion model with asymmetric potentials was developed 
to explain the bidirectional motions of kinesins. Based on this model, it was concluded that 
the hydrolysis of single ATP molecule is coupled to either forward or backward steps of 
kinesins. 

Although the theoretical picture presented by Nishiyama et al., 2002, that both forward 
and backward steps of kinesins are created by the same mechanochemical transduction 
mechanism, seems to be able to describe several features of the kinesin motility, there are 
serious fundamental problems with this view. It contradicts the current biochemical view 
of this process and earlier studies (Schnitzer and Block, 1997; Hua et al., 1997; Coy et al., 
1999) that show a tight coupling, i.e., one ATP molecule is hydrolyzed per each forward 
8-nm step. Note, however, that these earlier investigations mainly neglected the backward 
steps in their statistical analysis. In addition, the asymmetric potential used in the biased 
Brownian motion model breaks the periodic symmetry of the system, and it violates the 
principle of microscopic reversibility since the backward processes are not taken into account. 
Furthermore, this model cannot predict analytically the fraction of the forward and backward 
steps separately, and it also fails to account for irreversible detachments of kinesin molecules 
from the microtubules, which are observed in experiments. Clearly, a better quantitative 
theoretical description, which does not violate the basic physical and chemical principles, is 
needed in order to satisfactorily understand the mechanochemical coupling in kinesins. The 
aim of this article is to discuss in detail such a theoretical approach. 

We present a theoretical analysis of mechanochemical coupling and dynamics of ki- 
nesin molecules which utilizes the first-passage time processes (van Kampen, 1997) in pe- 
riodic discrete-state stochastic models. This is an extension of the recently developed ap- 
proach (Qian, 1997; Kolomeisky and Widom, 1998; Fisher and Kolomeisky, 1999a; Fisher 
and Kolomeisky, 1999b; Kolomeisky and Fisher, 2000a; Kolomeisky and Fisher, 2000b; 
Kolomeisky, 2001), which has been used successfully to analyze in detail the dynamics of sin- 
gle conventional kinesin molecules (Fisher and Kolomeisky, 2001) and myosin- V (Kolomeisky 
and Fisher, 2003). We argue that the experimental observations by Nishiyama et al., 2002, 
can be described by the simplest (A^ = 2)-state model with irreversible detachments, in 



which ATP hydrolysis is tightly coupled only to the forward steps of motor proteins. 

Theoretical Approach 
Chemical Kinetic Models 

Our approach is based on using multi-state discrete stochastic, or chemical kinetic, mod- 
els. The main assumption of the simplest periodic sequential chemical kinetic model, which 
is shown in Fig. la, is that a motor protein molecule is viewed as a particle that moves 
along a periodic linear track from one binding site to the next one through the sequence of 
N biochemical conformations. The particle in state j can jump forward to state j ' + 1 with 
the rate Uj, or it can slide one step backward to the site j — 1 with the rate Wj. After moving 
A^ sites forward the motor protein comes to the same biochemical state but shifted by a step 
size distance d. For kinesins this distance is 8.2 nm, and it is equal to the size of a tubulin 
subunit in microtubules (Howard, 2001). The states j = IN {I = 0, ±1, ±2, ■ ■ ■ ) represent 
the biochemical conformations where the motor protein molecule is tightly bound to the 
track, i.e., to the microtubule in case of kinesins, and without the ATP fuel molecule. ATP 
binding corresponds to the transitions from states j = IN to j = 1 + IN, while other forward 
transitions describe the ATP hydrolysis and subsequent release of hydrolysis products. It 
is important to note that, although the motor protein moves preferentially in one direction, 
the reverse transitions cannot be ignored in any reasonable model of motor protein motility, 
and the backward steps are frequently observed experimentally (Schnitzer and Block, 1997; 
Nishiyama et al., 2002). 

In the periodic sequential multi-state stochastic model the dynamics of the motor protein 
can be viewed as the motion of the particle on a periodic one-dimensional lattice (with a 
period A^). This observation allows one to derive an explicit analytical expressions for the 
mean velocity V{{uj,Wj}), 

V = hm ^, (1) 

in terms of transition rates {uj,Wj} for any value of A^ (Fisher and Kolomeisky, 1999a; 
Fisher and Kolomeisky, 1999b). Here, x{t) measures the position of the single molecule on 



the linear track. Specifically, the mean velocity is given by (Kolomeisky and Fisher, 2000a) 
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Note also the periodicity of transition rates, i.e., Uj±n- = Uj and Wj±n- = Wj. 

Similar arguments can also be applied to obtain closed-form exact analytic formulae for 
the dispersion D{{uj,Wj}) (or effective diffusion constant) of the motion, which is defined 
as follows, 

D = l \im ^[{x\t)) - {x{t)n (5) 

The simultaneous knowledge of both the velocity V and the dispersion D determines the 
bounds on rate-limiting biochemical transitions and thus provides a valuable information 
about the mechanism of motor proteins motility (Visscher et al., 1999; Kolomeisky and 
Fisher, 2000a; Fisher and Kolomeisky, 2001). 

One of the advantages of using chemical kinetic models to describe the processivity of 
motor proteins is the ability to easily incorporate the effect of external force F (Fisher and 
Kolomeisky, 1999, 2001). This can be done with the introduction of load-distribution factors, 
9j^ and 9j ( for j = 0, 1, ■ ■ ■ , iV — 1), that modify the transition rates in the following way, 

Uj => Uj{F) =u^jexp{-e+Fd/kBT), 

Wj => Wj{F) =w°jexp{e-Fd/kBT). (6) 

This is the consequence of the fact that the external load F modifies the activation barriers 
for forward and backward transitions, and the load-distribution factors refiect how they 
changed. It is also reasonable to assume that 

N-l 
j=0 

since the motor protein, making a step d against an external force F and going through N 
intermediate steps, produces a work equal to Fd. A force at which the motor protein stops 
moving is called a stall force. 



First-Passage Time Processes 

In many single-molecule experiments on motor proteins the fractions of forward and back- 
ward steps and dwell times between the consecutive events are measured (Nishiyama et al., 
2002; Asbury et al., 2003; Mehta et al. 1999). In terms of chemical kinetic models discussed 
above, these experimental quantities can be associated with the so-called splitting proba- 
bilities and conditional mean first-passage times, correspondingly. First-passage processes 
for sequential multi-state stochastic models are well studied (van Kampen, 1997; Pury and 
Caceres, 2003), and thus the available results can be easily adopted for the description of 
motor proteins dynamics. 

Consider a motor protein particle in state j, as shown in Fig. la. Recall that the 
sites — A^, and A^ correspond to the binding sites for motor proteins. Now let us define 
'^j,n{{uj, Wj}) as the probability that the particle starting from state j will reach the site A^, 
before backtracking to the previous binding site — A^. Similarly, we can define tt^ _Ar({-Uj, Wj}) 
as the probability for the particle to advance to state — A^ for the first time before reaching 
the forward binding site A^. These quantities are called the splitting probabilities (van 
Kampen, 1997). We are mainly interested in the case of j = 0, since the probabilities vro,±Ar 
give us the forward and backward fractions of stepping for the motor protein particle. The 
explicit expressions for splitting probabilities are known (van Kampen, 1997), and for the 
periodic A^-state stochastic models we obtain a simple relation, 

7ro,7v = 1 - vro,_7V = ^v-i. T^- (8) 

In a similar fashion, we can define the conditional mean first-passage times 
'Tj,±n{{uj,Wj}), that represent the average time the particle spends before advancing for- 
ward or backward to sites ±A^, correspondingly. Then it is easy to conclude that the dwell 
times for the forward steps of motor proteins correspond to To^n, while the dwell times for 
the backward steps are given by Tq^-at. The explicit expressions for the dwell times within 
the periodic A^-state chemical kinetic model can be derived from more general equations 
that are not restricted by periodicity conditions (see Pury and Caceres, 2003), yielding 

TTo.Af T^O-N /„v 

TO,N = , Tq-N = , (yj 

Ueff Weff 

where the effective transition rates -Ue// and We// are defined in Eqs. (3) and (4). Then 
applying the Eqs. (3) and taking into account the relations for the forward and backward 



fractions [see Eq. (8)], we conclude that 

ro,N = To-N- (10) 

This is a very important resuh because it indicates that the dwell times for the forward and 
backward steps are always equal to each other for any set of transition rates, although the 
probabilities of these steps may differ significantly. It is also important to note that periodic 
conditions in the system are crucial for this conclusion. 

Effect of Detachments 

Motor proteins do not always stay binded to the linear track, they can dissociate and 
diffuse away. For kinesins moving along the microtubules the effectively irreversible detach- 
ments have been observed experimentally (Schnitzer et al., 2000; Nishiyama et al., 2002). 
Theoretically, the effect of detachments on the drift velocity, dispersion and the stall force 
has been investigated (Kolomeisky and Fisher, 2000a) using the extension of the simplest 
sequential multi-state stochastic model. However, to best of our knowledge, the problem of 
how the motor protein dissociations change the mean first-passage time processes, namely, 
the fractions and mean dwell times of forward and backward steps, has not been studied 
at all. Below we outline how this effect can be solved by mapping it into another sequen- 
tial multi-state stochastic model but without detachments, for which the results are already 
known. 

Consider a motor protein particle in state j as shown in Fig. lb. It can move forward 
(backward) with the rate Uj (wj), or it can dissociate irreversibly with the rate 6j. We again 
define ttj^n and ttj^^n as splitting probabilities of reaching for the first time the forward (at 
A^) or the backward (at — A^) binding site. In addition, we introduce a new function iij^s 
as a probability for the motor protein, which starts at the site j, to detach before reaching 
the forward or the backward binding states. These probabilities are related through the 
normalization condition, 

TTj^N + TTj _Af + TTj^s = 1- (H) 

Now we may recall that the particle at the site j has to jump to the site j + 1 or j — 1, or 
it will detach. These jumps have the probabilities Uj/(uj + Wj + 6j), Wj/{uj + Wj + 6j) and 
^j/i'^j + "^i + '^i)' correspondingly. Then the expression for the forward splitting probability 



is given by (van Kampen, 1997) 

for any — A^ < j < N, and with the obvious choice of boundary conditions, 

7^N,N = 1, TTat-AT = 0. (13) 

Similar equations can be derived for the backward sphtting probabihties ttj^-n- 
The Eq. (12) can be easily rewritten as a difference equation, i.e., 

UjTTj^l^N + WjTCj^i^N - (Uj + Wj + Sj)nj^N = 0. (14) 

Assume that the solution of this equation can be presented as 

^j,N = 0i7r*7v, (15) 

where the function vr*^ is the splitting forward probability for a new system without de- 
tachments, and the auxiliary function (pj is yet to be determined. Substituting Eq. (15) into 
the Eq. (14) we obtain, 

Uj<f)j+lT^*j+l,N + Wj<Pj-l^j~l,N - (Uj + Wj + 5j)(l)j7ll^ = 0. (16) 

If we define new rates for the stepping process without detachments as 

u*j = Uj(j)j+i, Wj = Wj(f)j_i, (17) 

and also require that 

u* + w* = Uj(j)j+i + Wj(j)j^i = {uj + Wj + 5j)(j)j, (18) 

then the Eq. (16) is easily transformed into 

u*jT^*j+i,N + K^i-i,JV - («i + ^])t^In = 0, (19) 

with the boundary conditions tt^jv ~ ^ ^"^^ '^-n n ~ 0- These boundary conditions also 
mean that 0_Ar = (p^ = 1. Examining Eq. (19), one can observe that this is the expression to 
determine the forward splitting probability of the sequential multi-state stochastic process 
with rates {u*,Wj} but without detachments, for which the solutions are available (van 
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Kainpen, 1997). It leads to the explicit equation for the forward splitting probability. Similar 
arguments can be developed for the backward splitting probabilities. 

Our analysis relies on the ability to compute the functions (pj, which can be accomplished 
by utilizing the Eq. (18). However, it is more convenient to look at (j)j as elements of the left 
eigenvector of a (2A^ + l)x(2A^ + 1) matrix M, for which the non-zero elements are given 
by 

{uj + Wj + 5j), for i = j; 

Wi, for i = j — 1; (20) 

I Uj, for i = j + 1; 

with -A^ < i,j < N. 

The effect of detachments for conditional mean first-passage times can be investigated 
in a similar way. Here we again define r^ ^v {tj,-n) as the mean time to reach the forward 
(backward) binding state N {—N) for the first time. In addition, we define Tj^s as a mean 
first-passage time for the motor protein particle to dissociate from the molecular track before 
reaching the forward or backward binding sites ±A^. The mean first-passage times can be 
found by solving the backward Master Equation (see Pury and Caceres, 2003), 

UjTj+i^N + U!jTj_i^N - {Uj + Wj + 5j^N)rj,N = -1, (21) 

with the boundary conditions t±n^n = 0. Again, looking for the solution in the form 
'Tj,N = 'Pj'^j^Ni ^'^'i using Eqs. (17) and (18), we obtain the following expression, 

^j'^j+l,N + Wj^j~l,N - iUj + Wj)Tj^N = -1, (22) 

that determines the forward mean first-passage time for the system without detachments. 
Because the exact solutions for this case are available (van Kampen, 1997; Pury and Caceres, 
2003) the expressions for the mean first-passage times for the system with detachments can 
be easily obtained. 

The general equations for splitting probabilities and mean first-passage times are quite 
complex, and we present in the next subsection the expressions only for simple cases A^ = 1 
and N = 2. However, it can be shown that for any A^ the calculations of the mean dwell 
times to move forward, backward or to dissociate leads to the following important relation, 

T0,N = TQ-N = To,<5- (23) 

This is one of the main results of our theoretical analysis. 
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Results for A^ = 1 and N = 2 models 

To illustrate our method, let us consider two simple cases, A^ = 1 and N = 2 periodic 
sequential stochastic models with detachments. When the period of the system is A^ = 1, 
the auxiliary function 0o can be easily calculated, 

u + w ,^,, 

00 = ^-— T, 24 

U + W + 

and also recall that 0_i = 0i = 1. This leads to the simple relations for the splitting 
probabilities, 

T^o,! = u/{u + w + 6), TTo-i = w/{u + w + 6), tiq^s = S / {u + w + 6); (25) 

and for the mean first-passage times, 

ro,i = ro-i = ro,5 = 1/(m + w + (5). (26) 

For N = 2 case, the calculations become more tedious. The results for the functions 0_i, 
00, and 01 are given by 

UqUJ - UqwI + Wi{uq + Wq + 5q){ui + Wi + 5i) 



[{uq + Wo + 5q) {ui + wi + (5i) - (uqWi + UiWq)] (mi + Wi + 5i) ' 

, _ UqUi + WqWi 

[(Mo + Wq + 5o){Ui +Wi+ 5i) - {uqWi + UiWo)] ' 



(27) 



(28) 



(29) 



[{uq + Wo + 5o)(mi +wi + 6i) - (mo^i + uiWq)] (mi + Wi + 6i) ' 
Then, after lengthy but straightforward calculations, it can be shown that the splitting 
probabilities are 

[uqUi + wowi + 6q6i + 6q{ui + wi) + 6i{uq + wo)] ' 



WqWi 

[uqUi + WoWi + dodi + do(Mi + wi) + di[uo + Wo)J 
and -kq^s = 1 — 7ro,2 — t^q,-2- Similar calculations for the mean first-passage times yield 

_ _ _ Uq + Ui+Wq + Wi+5i 

To,2 — ^0,-2 — tqs — -, ■ -^—r- — -^"7 \ N , r / : rr- l<J^j 

\uqUi + woWi + dodi + do[ui + wi) + di[uo + wojj 
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Analysis of Kinesin Data 

Structural, biochemical and kinetic data on kinesins suggest that the protein molecule 
goes through at least four intermediate states (Lodish et al., 1995, Bray, 2001). However, 
recent study of kinesin dynamics using (A^ = 2)-state chemical kinetic model, which takes 
into account the irreversible detachments, provides a very reasonable description of some 
aspects of mechanochemical coupling in this system (Fisher and Kolomeisky, 2001). Thus, 
in order to analyze the experimental data of Nishiyama et al., 2002, we adopt the simplest 
model which includes only two states. The states j = ■ ■ ■ , —2, 0, 2, ■ ■ ■ would correspond to 
the kinesin with both molecular heads tightly bound to the microtubule and without ATP 
molecule. The states j = ■ ■ ■ , — 1, 1, ■ ■ ■ label all other kinesin conformations after ATP 
binding and subsequent hydrolysis and release of its products. 

It now follows that the forward ATP-binding transition depends linearly on ATP concen- 
tration, -Uq = A;q[ATP], where the superscript indicates the case of zero load: see also Eq. 
(6). At the same time the next forward rate Ui and the backward rate Wi do not depend on 
ATP concentration, while they may change under the effect of external forces. 

The final backward rate Wq might, in principle, depend on concentrations of ADP and 
inorganic phosphate that both are the products of ATP hydrolysis. However, most current 
experiments on kinesins utilize ATP regeneration system (Svoboda et al., 1994; Visscher et 
al, 1999; Schnitzer et al., 2000; Nishiyama, 2002), in which there are no independent control 
of [ADP] and [Pj]. As a result, we adopt a phenomeno logical description of this backward 
transition, namely, 

w', = K[ATP]/il + [ATF]/coY/', (33) 

where the parameter cq effectively describes the ATP regeneration process. This approach 
has been used successfully to describe the mechanochemical transitions in kinesin and 
myosin-V (Fisher and Kolomeisky, 2001; Kolomeisky and Fisher, 2003). Note, however, 
that the specific description of ATP regeneration process has a minimal effect in fitting of 
experimental results. 

After systematically exploring the multi-dimensional space of parameters and using Eqs. 
(30)-(33) the fractions of forward and backward steps and mean dwell times between the 
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consecutive steps of kinesins can be well described by the following rate constants 

k^ ~ 5.1fiM-^s-\ k'o ~ 2.8fiM'^s-\ Cq ~ 1.7/iM, w\ - 5.5 ■ 10~^s~\ 

ul ~ 121s~\ 5^ ~ l.ls'\ 5? ~ 1.6 • 10-V\ (34) 

and load-distribution parameters 

e+ ~ 0.0, e^ ^ 0.391, ^+ ~ 0.086, 
e^ ~ 0.523, el ~ 0.047, ^^ ~ 0.466. (35) 

The results of the fitting of experimental observations are given in Figs. 2 and 3. Note that 
the values for the parameters reported here are in a good agreement with the other indepen- 
dent investigation of kinesin motility (Fisher and Kolomeisky, 2001), where the multi-state 
periodic stochastic models have been used to analyze the single-molecule experimental mea- 
surements of velocities, stall forces and dispersions (Visscher et al., 1999). 

Discussion 

Our theoretical analysis provides explicit expressions for the fractions of forward and 
backward steps and dissociations, and for the mean dwell times between consecutive steps 
of motor proteins. This allows us to investigate the problem of mechanochemical coupling 
between the motion of kinesins and ATP hydrolysis. Our main conclusion is that the mean 
dwell times to move forward, backward or irreversible detach are equal to each other in- 
dependently of ATP concentration or external force. It means that the picture of tight 
coupling between ATP hydrolysis and forward steps of kinesins does not contradict the ex- 
perimental findings of Nishiyama et al., 2002. Moreover, the proposed bidirectional biased 
model (Nishiyama et al., 2002), which assumes that a hydrolysis of a single ATP molecule is 
coupled to either forward or the backward movement, is basically incorrect since it violates 
the principle of microscopic reversibility and breaks the symmetry of the system. 

Our theoretical results could also be understood in the following way. The mean dwell 
times between movements measured in single-molecule experiments are actually correspond 
to the mean lifetimes of states when the motor protein binds strongly to the linear track. 
Then these lifetimes should be independent of what direction the motor protein will go in 
the next step, although the probability of these steps might be rather different. 
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The analysis of mean dwell times at different external forces, as shown in Fig. 3, suggests 
that there is a maximum at high loads. This maximum is close but not exactly at the stall 
force. When [ATP] =10 ^M the maximum can be found at F ~ 6.6 pN, while the stall 
force is approximately equal to 6.8 pN. At high ATP (1 mM) the position of maximum is 
shifted to 7.7 pN, with the calculated stall force Fs — 9.2 pN. This can be understood in 
the following way. The external load decreases the forward transition rates, while slowing 
down the backward transitions. These two tendencies have an opposite effect on mean dwell 
times, and it leads to the observation of maximum at some specific value of external force. 

Because our method provides exact expressions for biophysical parameters, we are able 
to study the effects of ATP concentration and external forces on these parameters, and we 
can make a qualitative predictions that can be checked experimentally. 

First, we investigate how mean dwell times depend on [ATP] at different external loads. 
As shown in Fig. 4, the larger the external force, the larger is the mean dwell time. However, 
at constant force, the mean dwell time decreasing with increase in concentration of ATP. 
This is in agreement with intuitive expectations since at large [ATP] the binding process 
is faster. At the same time the external force slows down the binding and other forward 
processes more than it accelerates the backward transitions. These observations are also 
consistent with theoretical investigation of processivity of motor proteins using thermal 
ratchet approach (Parmeggiani et al., 2001). 

The dependence of the fractions of different movements on ATP concentration at different 
external loads is presented in Fig. 5. The increase in [ATP] increases the probability of the 
forward steps, while making the fractions of backward steps and detachments negligible. 
Finally, the predictions for the force-velocity based on the fitted parameters are given in 
Fig. 6. These predictions are generally agree with the values of drift velocities and stall 
forces obtained in other single-molecule experiments on kinesins (Visscher et al., 1999). 
However, the shapes for force-velocity curves are differ for ATP saturating conditions. 

Conclusions 

In summary, we have presented a theoretical study of mechanochemical coupling in ki- 
nesins. The analysis of multi-state stochastic models of motility using the method of first- 
passage times allowed us to obtain the explicit formulae for fractions of steps in different 
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directions, and for the mean dwell times between the steps, including the irreversible de- 
tachments. The experimental data on kinesins can be well described by this approach. Our 
analysis is consistent with the current theoretical view of tight coupling between catalytic 
cycles and mechanical steps for kinesins, i.e., one ATP molecule is hydrolyzed per each 
forward step, and the rare backward steps correspond to ATP production. Although our 
theoretical approach seems to provide a reasonable and convenient framework for investi- 
gating the mechanochemical coupling in different motor proteins, further experiments are 
needed in order to validate our theoretical picture. 
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FIGURE LEGENDS 

FIGURE 1. a) General schematic view of periodic multi-state stochastic models. A 
motor protein particle in state j can make a forward transition at rate Uj, or it can 
undertake a backward transition at the rate Wj. The states j = ■ ■ ■ ,—N,0,N,- ■ ■ 
correspond to the strongly bound states, b) General scheme of periodic multi-state stochas- 
tic models with irreversible detachments. The particle in state j can dissociate with a rate 6j. 

FIGURE 2. Probabilities, or fractions, of forward steps (circles), backward steps 

(triangles) and detachments (squares) as a function of the external force at a) [ATP]=lmM; 
b) [ATP] = 10/iM. 

FIGURE 3. Dwell times between the adjacent movements of the kinesin molecule 

as a function of external force. The filled symbols correspond to experimental measure- 
ments at [ATP] = 10/xM, while open symbols describe the experiments at [ATP]=lmM. The 
circles mark the experimental measurements for dwell times before the forward steps; the 
triangles correspond to experimental dwell times before the backward steps; and squares 
describe the dwell times before detachments. 

FIGURE 4. Predictions for the dwell times as a function of [ATP] at low {F = 1 

pN) and high external load {F = 5 pN). 

FIGURE 5. Predictions for the variation of the fractions of forward steps, back- 

ward steps and detachments at a) F = 1 pN; and h) F = 5 pN. 

FIGURE 6. Predictions for the force- velocity curves at different [ATP]. 
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